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摘 要 : 三 江 源 黄河 源 园区 作为 黄河 的 发 源 地 , 它 的 生态 系统 对 整个 黄河 流域 至 关 重 要 。 针 对 研究 
区 2000 一 2020 年 存在 的 土地 退化 问题 ,选取 土壤 侵蚀 模 数 、 沙 化 特征 指数 与 土壤 湿度 指数 ,采用 空 
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托 普 利 茨 逆 协 方差 聚 类 法 开展 了 3 个 参数 的 聚 类 分 析 , 从 而 进行 黄河 源 园区 土地 退化 时 空 演化 
分 析 , 并 采用 极致 梯度 提升 法 进行 时 空 演化 影响 因 


素 分 析 。 结 果 表 明 :(1) 土地 退化 现象 最 严重 的 


地 区 是 中 部 偏 南 地 区 与 东北 地 区 ,其 次 为 偏 北 地 区 。(2) 东北 地 区 存在 水 土 流失 与 土地 沙化 2 种 问 
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整体 好 转 , 偏 北 地 区 有 土地 沙化 的 风险 。 该 研究 对 


势 推进 黄河 源 园区 生态 保护 具有 重要 意义 。 


题 ,中 部 偏 南 地 区 以 严重 的 水 土 流 失 为 主 , 偏 北 地 区 主要 表现 为 土地 沙化 问题,(3) 近 20 a 来 ,土地 
退化 指数 呈现 波动 性 下 降 , 土 地 退化 整体 趋 于 好 转 , 在 偏 南 地 区 表现 得 更 明显 ,(4) 水 土 流失 问题 
| 掌握 黄河 源 园区 土地 退化 的 时 空 分 布 及 发 展 赵 


X 键 词 : 土地 退化 ; ER; 空间 托 普 利 茨 逆 协 方差 聚 类 ; 极致 梯度 提升 
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土地 退化 被 定义 为 土地 生产 能 力 的 暂时 或 永 
入 衰 退 ”, 主要 分 为 土壤 物质 的 迁移 和 土壤 内 部 化 
学 性 质 的 变化 ”。 由 于 其 对 环境 .生物 与 农业 等 方 
面 带 来 的 严重 破坏 作用 ,得 到 国内 外 相关 部 门 广泛 重 
视 。 然 而 ,对 土地 退化 的 评 佑 与 监测 尚未 达到 一 个 共 
YA ,特别 是 面向 区 域 及 以 上 空间 尺度 的 研究 六 。 

土地 退化 问题 研究 主要 以 遥感 手段 为 主 , 通 过 
遥感 反 演 建立 土地 退化 指数 "的 方法 进行 土地 退 
化 评价 。 土 壤 “ 土地 利用 /土地 覆 被 ”植被 、 
地 表 反 照 率 ”” 气候 ”等 因子 均 被 前 人 应 用 于 建立 
土地 退化 指数 。 近 年 来 ,为 了 得 到 更 准确 的 评估 、 
监测 与 预测 结果 ,土地 退化 研究 更 加 趋向 于 采用 多 
源 时 空 数据 ,应 用 机 顺 / 次 度 学 习 算 法 充分 挖掘 多 源 
时 空 数据 中 隐 含 的 信息 。Yue 等 ?5 采用 反 向 传播 神 
经 网 络 (Back propagation-artificial neural network , BP- 
ANN) 得 到 了 鄂尔多斯 高 原 的 土地 退化 监测 结 
许 明 杰 等 以 一 种 加 权 自 组 织 映 射 神经 网 络 (Self or- 
ganizing map,SOM) 对 黄冈 市 土地 生态 进行 了 评估 ; 
Grinand 等 ”采用 随机 森林 算法 进行 了 土地 退化 建 模 
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与 预测 。 研 究 结果 表明 ,机 带 / 深 度 学 习 算 法 在 土地 
退化 分 类 与 预测 方面 均 取 得 了 较 好 的 效果 。 但 
是 ,综合 时 空 数据 进行 土地 退化 时 空 演化 分 析 的 研 
究 仍 有 不 足 , 结 合 地 理 空 间 人 工 智 能 (Ceospatial ar- 
tificial intelligence , Geo AI) 进行 后 续 研 究 则 更 为 匮 

三 江 源 黄 河源 园区 作为 黄河 的 源头 ,对 整个 黄 
河流 域 至 关 重 要 ,但 由 于 水 土 流失 与 土地 沙化 等 问 
题 的 存在 导致 生态 发 展 受 限 ”" ,国家 对 此 也 越 来 
越 重视 。 特 别 是 近 20 a 以 来 ,为 了 保护 这 片 敏感 的 
地 区 《青海 三 江 源 自然 保护 区 生态 保护 和 建设 总 
体 规划 》 等 政策 相继 出 台 。 为 了 响应 国家 的 号 召 ， 
本 文 旨 在 对 当地 的 土地 退化 治理 提供 理论 依据 。 根 
据 研究 区 现存 的 水 土 流失 与 土地 沙漠 化 问题 “” ， 
以 2000.2005.2010.2015 年 和 2020 年 5 期 的 土壤 侵 
蚀 模 数 .沙化 特征 指数 与 土壤 湿度 指数 为 评价 指 
标 ,应 用 空间 托 普 利 茨 逆 协 方差 聚 类 (Spatial To- 
eplitz inverse covariance-based clustering, STICC ) 进 


行 土地 退化 空间 分 布 模式 分 析 ,并 结合 土地 退化 指 
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数 分 析 其 时 空 演化 过 程 。 在 此 基础 上 ,采用 极致 梯 
度 提升 (Extreme gradient boosting, XGBoost) 算 法 
进行 土地 退化 时 空 演化 驱动 机 制 探索 。 


1 研究 区 概况 


= YT YG BE Yay VE BK (97° 01' 20" ~99° 14'57"E, 
33°55'05"~35°28'15"N ) 位 于 青海 省 果 洛 藏族 自治 州 
玛 多 县 境内 "5( 图 1) ,地 处 青藏 高 原 ,海拔 4125~ 
5257 m, 中 部 偏 南 与 东北 地 区 以 山地 地 貌 为 主 , 西 
部 与 北部 地 区 以 湖泊 湿地 为 主 ,拥有 扎 陵 湖 、 鄂 陵 
湖 2 个 最 大 的 天 然 湖泊 。 植 被 主要 分 布 在 研究 区 东 
南 、 西 南 地 区 ,北部 相对 稀 琉 。 气 温 .降雨 年 度 差异 
不 大 ,大 致 呈现 由 南 向 北 逐 渐 降低 的 趋势 。 年 降雨 
量 100~200 mm ,集中 在 6 一 10 月 ,年 平均 气温 维持 
在 1? 左 右 。 但 是 随 着 全 球 气 候 变 化 高原 生态 环境 
暧 干 化 趋势 加 剧 ,再 加 上 人 类 不 合理 活动 的 双重 作 
用 ,三 江 源 黄河 源 园 区 生态 环境 日 益 恶 化 水土 流 
失 加 剧 、 土 地 沙漠 化 面积 扩大 1。 


图 例 
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图 1 三 江 源 黄 河源 园区 地 形 图 
Fig. 1 Topographic map of Sanjiangyuan Yellow River 
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2 数据 与 方法 


2.1 数据 来 源 

数字 高 程 模型 (Digital elevation model, DEM ) 数 
据 : 该 数据 来 源 于 美国 国家 航天 局 地 球 科学 数据 系 
统 (http://Earthdata.nasa.gov/)2019 年 发 布 的 DEM V3 
数据 ,其 空间 分 辩 率 为 30m ,垂直 精度 为 20 m。 

降雨 数据 :降雨 数据 来 源 于 ERA5 数据 集 ” ,由 
于 2020 年 的 降雨 数据 不 全 ,本 文 基于 该 数据 集 选 取 
了 2000、2005 .2010、2015 年 和 2019 年 共 5 a 的 逐日 
降雨 量 数据 。 

土壤 数据 :土壤 数据 来 自 于 中 国 科 学 院 资源 环 


境 科 学 数据 中 心 (http://www.resdc.cn/) 的 1 km 分 辩 
率 的 土壤 质地 类 型 数据 ,该 数据 集 采 集 于 全 国 第 二 
次 土壤 普查 期 间 ;以 及 国家 青藏 高 原 科学 数据 中 
心 (http:/data.tpde.ac.cn/) 的 中 国土 壤 有 机 质 食量 
数据 “ ,其 元 数据 更 新 于 2021 年 。 

植被 覆盖 度数 据 : 植 被 覆盖 度数 据 来 源 于 
MOD13A1 的 NDVI 数 据 集 , 由 于 7 一 9 月 中 旬 植 被 旺 
盛 , 因 此 分 别 统计 2000、2005、2010、2015 年 和 2020 
年 Sa 的 NDVI 均 值 。 并 在 Google Earth Engine 平台 
上 ,通过 像 元 二 分 法 进行 处 理 后 得 植被 覆盖 度 

遥感 数据 :利用 Landsat 系列 卫星 的 蓝 (BLUE)、 
短波 红外 1(SWIR1) 以 及 短波 红外 2(SWIR2) 提 取 
沙化 特征 指数 。 为 了 避免 植被 与 南 雪 的 干扰 ,在 影 
像 选取 时 以 10 一 11 月 为 时 间 区 间 ,以 云 量 低 于 15% 
为 宜 。 由 于 2015 年 和 2017 年 的 Landsat 7 影像 存在 
损坏 ,2015 Æ AY Landsat 5 和 2020 年 的 Landsat 8 $% 
像 在 该 时 间 段 内 云 量 较 高 ,最终 选择 了 2000 年 的 
Landsat 7.2004 年 和 2010 年 的 Landsat 5.2015 年 和 
2019 年 的 Landsat 8 ALMA, BSP RS 
提出 的 土壤 湿度 指数 提取 方法 ,选取 MOD09A1 E. 
星 中 的 第 1~7 波 段 提取 土壤 湿度 指数 。 由 于 北方 每 
年 7 一 9 月 降雨 较 充分 .土壤 湿度 较 高 , 且 2000 年 的 
卫星 影像 存在 异常 ,分 别提 取 2001 , 2005, 2010, 
2015 年 和 2020 年 7 一 9 月 中 名 土壤 湿度 指数 均值 作 
为 研究 数据 。 
2.2 研究 方法 
2.2.1 土地 退化 评价 指标 ”参考 前 人 研究 成 果 *"*"， 
选用 土壤 侵蚀 模 数 表征 水 土 流失 问题 ,以 沙化 特征 
指数 与 土壤 湿度 指数 探究 土地 沙化 问题 。 共 线性 
检验 得 任何 2 个 因子 之 间 的 方差 膨胀 系数 (Vari- 
ance inflation factor, VIF ) 均 处 于 10 以 下 ,所 选 因 子 
之 间 不 存在 多 重 共 线 性 问题 。 

(1) 土壤 侵蚀 模 数 

自 20 世 纪 起 ,国内 外 研究 广泛 采用 土壤 侵蚀 模 
数 来 反映 水 土 流失 问题 。 在 区 域 尺度 下 ,多 采用 通 
用 土壤 流失 方程 (Universal soil loss equation , USLE) 
进行 水 土 流失 监测 ” ,具体 计算 公式 如 下 : 

A=RXKXLSXCxP (1) 

式 中 : 4 为 土壤 侵蚀 模 数 (1.km”…a');R 为 降雨 侵蚀 
力 因 子 (MJ.mnm.hnm2 .ha ;天 为 土壤 可 蚀 性 因子 
(thmz.hho2.MJT mm ;LS 为 坡 长 坡度 因子 ;C 
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为 植被 覆盖 因子 ;P 为 水 土 保持 工程 措施 因子 。 

在 实际 计算 过 程 中 ,可 根据 不 同 地 区 的 地 理 条 
件 进行 相应 修正 。 根 据 降雨 数据 ,R 值 计算 采 
用 音 文 波 等 ”提出 的 计算 方法 。 根 据 张 科 利 等 
针对 中 国 地 区 进行 了 修正 后 的 K 值 计算 方法 ,基于 
土壤 数据 进行 土壤 可 人 刨 性 计算 。 以 DEM 数据 为 准 ， 
在 LS 因子 计算 过 程 中 ,缓坡 坡度 因子 采用 MeCool ” 
的 方法 ,陡坡 坡度 因子 采用 Liu 等 ”提出 的 方法 , 结 
A Foster 等 中 提出 的 非 流 量 累 积 坡 长 计算 方法 与 
Wischmier 等 ”提出 的 坡 长 因子 计算 公式 提取 坡 长 
因子 。 基 于 植被 覆盖 度数 据 , C (ER RIE 
针对 中 国 地 区 提出 的 方法 。P 值 采用 章 影 等 “使 用 
的 方法 。 

(2) 沙化 特征 指数 


XX 
X ae = as ae X. (5) 
UH: Xo Xa 分别 为 单 因子 正 向 38 e H — 4e 
后 的 结果 ; X HERRE; X a 为 5 期 数据 中 的 最 大 
值 ; Xun 为 5 期 数据 中 的 最 小 值 。 
2.2.2 土地 退化 室 间 分 布 模式 

(1) 土地 退化 评价 单元 提取 

流域 作为 地 貌 的 最 小 单元 ,最 合适 作为 土地 退 
化 的 评价 单元 。 为 了 提取 到 合适 的 流域 ,需要 找到 
最 佳 汇 水 阔 值 汪 。 并 基于 该 阔 值 ,以 DEM 为 数据 
源 , 通 过 水 文 分 析 进 行 流域 提取 。 前 人 研究 表明 ， 
通过 对 不 同 汇 水 阔 值 下 的 河 网 密度 进行 均值 变 点 
分 析 , 求 取 变 点 后 ,可 得 到 最 佳 汇 水 阔 值 (图 2)。 分 
析 可 得 , 汇 水 国 值 为 300 的 点 为 变 点 ,以 此 将 整个 黄 


由 于 已 沙化 地 物 的 光谱 特征 与 其 他 地 物 存 在 
明显 差别 ,可 根据 Landsat 系列 卫星 的 BLUE、SWIR1 
及 SWIR2 波段 反 演 得 到 土地 沙化 程度 5 E 
Google Earth Engine 平 台 上 对 Landsat 系 列 数据 集 进 
行 去 云 .拼接 .裁剪 以 及 波段 运算 后 提取 沙化 特征 
指数 ,计算 公式 如 下 : 

SFI =(SWIR1 - BLUE)/(200 - SWIR1 - SWIR2) (2) 
式 中 :SFI 为 沙化 特征 指数 ;BLUE SWIR1 与 SWIR2 
分 别 对 应 着 Landsat 的 蓝 、 热 红外 1 及 热 红 外 2 波段 。 

(3) 土壤 湿度 指数 

土壤 湿度 指数 既 反 映 了 十 地 干旱 程度 ,又 反映 
了 土地 沙化 的 风险 55。MODIS 作为 具有 高 时 间 分 
辩 率 的 高 光谱 数据 ,在 反 演 土壤 湿度 时 起 到 了 较 好 
AY 3 32", Æ Google Earth Engine 平台 上 对 
MOD09A1 Zhe AE RF aM ae HTT EE PEE BE A 
计算 公式 如 下 : 

WET=0.0001 x(0.1147 X b1 +0.2489 x b2 + 


0.2408 x b3 + 0.3132 x b4 - 0.3122 x b5 - (3) 
0.6416 x b6 - 0.5087 x b7) 


st AO: WET X E HE i E FG BX 5 b1~b7 分 别 为 
MODO9AL 卫星 的 波段 1 到 波段 7。 

(4) 评价 指标 归 一 化 

由 于 所 选 土地 退化 影响 因子 存在 量 纲 不 一 致 
性 ,因此 需 对 其 进行 归 一 化 。 土 壤 侵蚀 模 数 与 沙化 
特征 指数 是 正身 指标 , 需 进 行 正 向 归 一 化 ;土壤 湿 
度 指 数 是 逆向 指标 , 需 进 行 逆向 归 一 化 。 具 体 归 一 
化 公式 如 下 : 


A = -i (4 ) 
min 


河源 园区 划分 为 44816 个 流域 ,并 将 此 作为 土地 退 
化 空间 分 布 模式 分 析 的 基础 。 


Ñ N, 
pe ~ 
Ray 
i ~ 
& Ss 
区 0 ` 
E Ss 
E Sie 
10 L 1 1 J 
0 500 1000 1500 2000 
汇 水 阔 值 /个 


图 2 Ae Ek Bete HER 


Fig. 2 Extraction of optimal catchment threshold 


(2) 土地 退化 空间 分 布 模式 

为 了 更 优 地 显示 土地 退化 空间 分 异 现象 ,开展 
了 SOM 神 经 网 络 .基于 密度 的 聚 类 方法 (Density- 
based spatial clustering of applications with noise ,DB- 
SCAN) AI STICC 聚 类 算法 对 比分 析 。STICC 聚 类 结 
果 的 空间 连续 性 强 于 基于 属性 的 SOM 神经 网 络 , 相 
似 属 性 对 象 的 聚集 优 于 基于 密度 的 DBSCAN。 可 
见 ,STICC 聚 类 方法 更 适合 进行 地 理 现 象 空间 分 布 
问题 的 研究 。 该 方法 采用 Delaunay = ff} Fill IE HR 
持 了 空间 上 的 连续 性 ,使 用 多 源 高 斯 分 布 的 逆 协 方 
差 矩 阵 定义 子 区 域 属 性 组 成 的 马尔 科 夫 随机 场 , 同 
时 考虑 空间 连续 性 与 重复 地 理 模 式 发 现 问题 ”。 

考虑 到 样本 数据 之 间 的 不 确定 性 与 随机 性 存 
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在 差异 ,采用 炉 值 法 ”对 正 向 逆向 归 一 化 后 的 土 
壤 侵 蚀 模 数 沙化 特征 指数 与 土壤 湿度 指数 分 别 确 
权 后 ( 表 1), 以 44816 个 流域 作为 聚 类 评价 单元 进行 
样本 值 的 提取 ,最 终 得 到 一 个 3x44816 的 矩阵 ,并 以 
此 进行 STICC 聚 类 。 为 了 避免 异常 值 对 聚 类 结果 的 
于 扰 ,首先 采用 局 部 Moran s IER EIT 
值 的 去 除 。 聚 类 时 ,将 分 类 个 数 统一 设置 为 3 ,将 逆 
协 方差 矩阵 稀 朴 性 水 平 调整 指数 (lambada) 统 一 设 
置 为 10” ,将 空间 一 致 性 调整 指数 (beta) 统 一 设置 为 
102 ,迭代 次 数 分别 设 置 为 5 .10、…、30。 为 了 得 到 
最 佳 迭 代 次 数 , 采 用 分 类 适 确 性 指数 (Davies-Bould- 
in index,DBI) 指 数 ” 对 聚 类 结果 进行 评价 ,DBI 指 数 
越 小 表明 聚 类 效果 越 好 。 结 合 DBI 指 数 得 到 ,2000、 
2005 , 2010, 2015 年 和 2020 年 的 迭代 次 数 分 别 为 
25 30.10.15 和 30 时 聚 类 结果 较为 适宜 。 


表 1 土地 退化 指数 权重 确定 


Tab. 1 Weight determination of land degradation indices 


指数 2000 年 2005 年 2010 年 20154F ” 2020 年 
土壤 侵蚀 模 数 0.768 0.504 0.606 0.887 0.870 
( 正 向 指标 ) 
沙化 特征 指数 0.162 0.454 0.329 0.025 0.044 
( 正 向 指标 ) 
土壤 湿度 指数 0.070 0.042 0.065 0.089 0.086 
(逆向 指标 ) 


| 


2.2.3 土地 退化 时 宝 演 化 探索 ”由 于 聚 类 结果 仅仅 
能 保证 将 相似 的 样本 分 为 一 类 ,无 法 得 到 该 类 所 代 
表 的 土地 退化 程度 高 低 , 也 无 法 在 时 空 尺 度 上 对 土 
地 退化 程度 进行 比较 。 为 了 探索 土地 退化 时 空 演 
化 过 程 ,还 需要 结合 土地 退化 指数 。 该 指数 分 布 在 
0~1 之 间 , 值 越 高 ,土地 退化 程度 越 高 。 整 合 5 期 土 
地 退化 指数 数据 ,采用 自然 间断 点 分 级 法 可 将 土地 
退化 程度 分 为 5 级 ,土地 退化 指数 处 于 0~0.098 为 微 
度 土 地 退化 .0.098~0.162 为 轻 度 土地 退化 .0.162~ 
0.217 为 中 度 土 地 退化 .0.217~0.290 为 微 重 土地 退 
化 、>0.290 为 重度 土地 退化 。 

LDI = a x USLE, „ie + b X SFL +cxWET (6) 
式 中 : LDI 为 土地 退化 指数 ; USLE ae . SEke WIE 
向 归 一 化 土壤 侵蚀 模 数 .沙化 特征 指数 ; WET 为 
道 向 归 一 化 土壤 湿度 指数 ;a、b 及 ce 分 别 代表 对 应 
权重 。 

综合 分 析 5 期 STICC 聚 类 结果 中 3 类 的 类 别 内 
土地 退化 指数 均值 ,可 将 3 类 进一步 定性 评价 为 高 、 


中 值 及 低 值 聚 类 区 。 在 此 基础 上 分 别 统 计 5 期 
STICC 聚 类 结果 中 每 一 类 内 的 土地 退化 指数 分 级 结 
果 面 积 占 比 ,根据 面积 占 比 随时 间 演 变 的 结果 可 进 
行 土地 退化 时 空 演化 过 程 探 索 。 

2.2.4 土地 退化 驱动 机 制 分 析 为 了 深入 探究 土地 
退化 时 空 演 化 驱动 机 制 ,本 文采 用 XGBoost 算 法 ， 
以 每 一 类 中 的 正 向 . 闭 向 归 一 化 土壤 侵蚀 模 数 、 沙 
化 特征 指数 与 土壤 湿度 指数 与 LDI 分 级 结果 为 数据 
源 , 对 不 同类 进行 特征 重要 性 评估 。 


3 结果 与 分 析 


3.1 土地 退化 空间 分 布 模式 分 析 及 结果 验证 

STICC 聚 类 将 整个 研究 区 划分 为 3 类 ,第 工 类 
分 布 面积 最 少 , 仅 占 10% 左 右 ,处 于 中 部 偏 南 与 小 
部 分 东北 地 区 ; 除 2005 年 外 ,第 下 类 占 整个 研究 区 
面积 大 约 50% ,大 体 分 布 在 偏 北 地 区 ;第 亚 类 除 
2005 年 面积 占 比 较 高 外 ,其 余年 份 占 整个 研究 区 面 
积 比 约 为 40%, 以 西南 、 东 南 地 区 分 布 最 广 (图 3)。 

第 [类 所 处 区 域 土地 退化 程度 相对 较 高 ,属于 
高 值 聚 类 区 ;第 开 类 为 中 值 聚 类 区 ;第 亚 类 表现 为 
低 值 聚 类 区 (图 3)。 结 合 2020 年 Landsat 8 真 彩色 
Fr BABE” 与 2020 年 的 聚 类 结果 进行 对 比 后 得 到 
实际 的 土地 退化 情况 与 真实 遥感 反 演 的 结果 基本 
一 致 。 示 例 3 和 5 表现 为 土地 退化 程度 轻 或 者 不 存 
在 土地 退化 的 一 类 ,植被 覆盖 度 高 ,整体 呈现 为 绿 
色 或 者 深 绿色 ;示例 1.2 和 6 基本 处 于 土地 退化 中 
值 聚 类 区 ,植被 覆盖 度 低 、 以 裸 土 分 布 为 主 ,地 表 呈 
现 为 较 亮 的 颜色 ,地 表 分 布 较为 平整 ,整体 色差 不 
大 ;示例 4 和 7 表现 为 土地 退化 程度 最 重 的 一 类 ,地 
表 破 碎 、 沟 蜜 纵横 ,土壤 颜色 偏 红 、 偏 黑 ( 图 4)。 所 
以 利用 STICC 聚 类 进行 土地 退化 评价 是 可 靠 的 。 
3.2 土地 退化 时 空 演 化 探索 分 析 

基于 空间 分 布 模式 分 析 结 果 ,分 别 统 计 5 期 
STICC 聚 类 结果 中 第 [ TKR BMRA AHIR 
化 指数 分 级 结果 (图 3、 表 2)。2000 年 土地 退化 程度 
以 轻 度 为 主 , 占 整个 研究 区 面积 的 92.68%。 第 焉 类 
占 整个 研究 区 面积 的 52.66% , 主要 分 布 在 西南 地 
区 ,在 东南 地 区 也 有 少量 分 布 ;第 开 类 分 布 在 偏 北 
地 区 ,占据 整个 研究 区 面积 的 36.56% ;第 工 类 主要 
分 布 在 中 部 偏 南 地 区 与 东北 偏 北 地 区 ,但 分 布 面积 
仅 占 整个 研究 区 面积 的 10.78% , 且 以 轻 .中 度 为 
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(a) 20004 (b) 200547: 


À 


(c) 20104 


>z 
>z 


~ 0 25km 
— | 


35H ALDL E 类 别 内 LD 均值 Q 5km a MLD ae 
mm 0.146 0.132 mm 0.126 mm 0.275 0.252 mm 0.238 mm 0.208 0.202 mm 0.197 
(d) 201545 (e) 20204F 
N 
A 图 例 
O 黄河 源 园区 边界 
水 域 及 其 边界 
STICC 到 类 结果 
me 1% 
g 第 I 类 
类 别 内 LDI 岁 什 号 sm KLD Eo mm $ 
mm 0.079 0.071 mm 0.066 mm 0.088 0.075 mm 0.069 
注 :LDI 为 土地 退化 指数 。 
图 3 土地 退化 空间 分 布 
Fig. 3 Spatial distribution of land degradation 
(g) 示例 7( 第 工 类 ) 


(a) 示例 1( 第 工 类 ) 


( aN Blo M38) 


例 3( 第 亚 类 ) 4) 示例 4( 第 I 28) (e) RBIS IS 


图 4 土地 退化 空间 分 布 结 果 对 比 


Fig. 4 Comparison of spatial distribution results of land degradation 


主 。2000 一 2005 年 ,土地 退化 程度 由 以 轻 度 为 主 转 
化 为 以 微 重 为 主 ,土地 退化 程度 加 深 , 第 陡 类 大 幅 
增长 ,第 工 类 与 第 工 类 则 存在 一 定 程 度 的 减少 。 这 
表明 土地 退化 程度 在 偏 西 南 地 区 增长 趋势 高 于 其 
他 地 区 ,第 工 类 中 的 中 部 偏 北 地 区 与 西南 地 区 之 间 


的 差异 减 小 ,中 部 偏 北 地 区 转化 为 第 五 类 ,第 五 类 
土地 退化 面积 增长 了 33.29%。 到 2010 年 时 ,土地 
退化 程度 由 以 微 重 为 主 好 转 至 以 中 度 为 主 , 89.28% 
的 区 域 处 于 中 度 土 地 退化 。 由 于 偏 南 地 区 的 土地 
退化 好 转 趋势 相对 较 快 ,第 亚 类 中 的 偏 北 与 偏 南 地 
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表 2 土地 退化 程度 时 空 演化 统计 


Tab. 2 Spatial and temporal evolution statistics of land degradation degree 


面积 占 比 /% 
2000 年 2005 年 2010 年 2015 年 2020 年 
第 [类 第 [I 类 PIX 第 I 类 第 [I 类 第 于 类 第 I 类 第 [类 第 亚 类 第 I 类 第 [类 BM 第 I 类 第 [类 HIX 
WUE 001 0.00 0.01 0.00 0.00 0.00 0.00 0.00 0.00 11.35 41.58 41.19 5.25 62.33 23.12 


轻 度 7.53 35.04 50.11 0.00 0.00 0.00 0.00 
中 度 2.39 1.41 1.82 0.12 0.14 1.94 14.67 
微 重 0.61 0.10 0.52 3.81 6.95 81.54 4.42 


重度 0.24 0.01 0.20 2.02 1.01 2.47 0.08 


区 之 间 的 差异 逐渐 增 大 , 偏 北 地 区 的 第 焉 类 向 第 
1 工 类 转化 ,第 开 类 面积 大 幅 增长 ,占据 整个 研究 
区 面积 的 56.93%。 同 时 , 偏 南 地 区 的 第 I 类 由 于 土 
地 退化 指数 的 降低 ,逐渐 与 偏 北 、 偏 东 地 区 的 一 部 
分 第 开 、 焉 类 融合 ,面积 增长 了 13.22% ,但 第 I 8+ 
地 退化 程度 好 转 至 以 中 度 为 主 。2010 一 2015 年 , 土 
地 退化 指数 继续 下 降 ,并 转化 为 以 微 度 为 主 ,处 于 
微 度 的 区 域 占据 整个 研究 区 的 94.12%。 土 地 退化 
好 转 趋势 由 南 向 北 扩展 ,由 于 中 部 略 偏 北 土地 退化 
现象 好 转 , 有 18.79% 的 中 部 偏 北 地 区 的 第 TS 
转化 为 第 亚 类 。2015 一 2020 年 ,土地 退化 程度 基本 
稳定 ,继续 以 微 度 为 主 , 但 第 下 类 面积 增加 ,5.71% 
的 第 工 类 与 17.54% 的 第 焉 类 转化 为 第 下 类 ,这 主要 
是 由 于 中 部 偏 西 的 第 亚 类 土地 退化 指数 存在 小 幅 
度 的 提高 .中 部 偏 东 的 第 I 类 土地 退化 指数 存在 小 
幅度 的 下 降 。 同 时 ,由 于 东北 地 区 的 土地 退化 指数 
存在 小 幅 上 涨 ,3.42% 的 区 域 土地 退化 程度 加 重 至 
微 度 以 上 。 
3.3 土地 退化 驱动 机 制 分 析 

本 文采 用 XCBoost 算 法 呈 ,对 5 期 数据 的 3 个 类 
别 分 别 进行 土地 退化 特征 因子 重要 性 分 析 ( 图 5) 发 
现 , 土 壤 侵 蚀 模 数 是 造成 土地 退化 程度 空间 分 异 的 
主要 原因 ,但 2005 年 的 第 五 类 、2015 年 和 2020 年 的 
第 工 类 土壤 湿度 指数 对 土地 退化 空间 分 异 的 影响 
力 有 所 提高 。 结 合 单 因子 分 布 情况 (图 6) 可 知 ， 
2000 一 2005 年 土地 退化 程度 加 重 , 水 土 流失 风险 提 
高 土地 沙化 风险 小 幅 降低 。 由 于 降雨 侵蚀 力 由 
19.92~276.79 MJ mm hm .ha 增长 至 76.39~ 
695.06 MJ:mm:hm*…h'…a', 且 由 中 部 偏 南 地 区 向 偏 
西南 地 区 偏 移 ,植被 覆盖 度 仅 在 中 部 偏 北 地 区 有 小 
幅 增长 。 因 此 ,土壤 侵蚀 模 数 在 中 部 偏 南 、 特 别 是 
偏 西南 地 区 的 土壤 侵蚀 模 数 大 幅 增 长 ,逐渐 趋 近 于 


0.00 0.00 2.34 1.52 1.49 2.34 3.67 1.73 
52.39 22.22 0.23 0.09 0.14 0.35 0.36 0.31 
4.49 1.63 0.03 0.00 0.05 0.27 0.06 0.11 
0.05 0.04 0.00 0.00 0.00 0.03 0.01 0.06 


土壤 湿度 指数 国 沙化 特征 指数 国 土壤 侵蚀 模 数 
IIM II 亚 INM II 了 ITM 


U 


影响 因子 重要 性 


0.0 l 1 | | | 
2000 2005 2010 2015 2020 


年 份 


图 5 土地 退化 影响 因子 重要 性 


Fig. 5 Importance of land degradation impact factors 


土地 退化 程度 较 高 的 偏 北 地 区 ,第 焉 类 分 布 面积 

幅 增 加 。 同 时 ,土壤 湿度 在 中 部 偏 东 南 地 区 有 所 增 
长 ,但 是 在 中 部 偏 北 地 区 变化 不 大 ,土壤 湿度 指数 
对 第 亚 类 的 土地 退化 程度 空间 分 异 的 影响 力 大 幅 
提高 。 偏 东南 地 区 的 土地 沙化 风险 降低 ,土地 退化 
程度 加 重 趋 势 相 对 较 缓 。2010 年 ,降雨 侵蚀 力 在 中 
部 偏 西 南 地 区 的 高 值 由 695.06 MJ .mm .hm .ha 
下 降 至 268.59 MJ.mm:hm*…h'…a', 偏 西北 地 区 的 
降雨 侵蚀 力 相 对 提高 ,植被 覆盖 度 在 中 部 偏 北 地 区 
轻微 下 降 。 因 此 , 偏 南 地 区 的 土地 退化 好 转 趋势 高 
于 北方 ,中 部 偏 西 南 地 区 与 偏 北 地 区 的 空间 分 异 拉 
大 ,分 异 成 2 类 ,这 也 表明 了 2005 年 中 部 偏 西 南 地 
区 的 土壤 侵蚀 模 数 提高 并 不 具备 长 期 发 展 趋 势 。 
同时 , 偏 北 地 区 沙化 特征 指数 下 降 .土壤 湿度 提高 ， 
土地 沙化 风险 降低 。2010 一 2015 年 ,中 部 地 区 的 植 
被 覆盖 度 显 车 提高 ,土壤 侵蚀 模 数 继续 保持 下 降 趋 
势 。 但 偏 北 地 区 的 土壤 湿度 下 降 .沙化 特征 指数 提 
高 ,土地 沙化 风险 加 大 ,因此 第 开 类 分 布 向 北 偏 
移 。 第 下 类 中 的 偏 东 地 区 土地 沙化 现象 更 严重 , 造 
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(a) 土壤 侵蚀 模 数 (2000 年 ) (b) 土壤 侵蚀 模 数 (2005 年 ) (c) 土壤 侵蚀 模 数 (2010 年 ) (dA 土壤 侵蚀 模 数 (2015 年 ) (e) 土壤 侵蚀 模 数 (2020 年 ) 


N z ~ 
人 =A 4 


be Ga 
L 9 50km 9 50km 9 50km 5 
土壤 侵蚀 模 数 土壤 侵蚀 模 数 土壤 侵蚀 模 数 土壤 侵蚀 模 数 土壤 侵蚀 模 数 
0.00 WM 0.49 0.00 E = 1.00 0.00 © = 0.28 0.00 F 0.25 0.00 f= = 0.36 


(A 沙化 特征 指数 (2000 年 ) (e) 沙化 特征 指数 (2005 年 ) (h) 沙化 特征 指数 (2010 年 ) ”GD 沙化 特征 指数 (2015 年 ) ”0) 沙化 特征 指数 (2020 年 ) 


Ce ee ee Oe oe 


Q {0m 9 50km 9 50km 9 50km 9 50km 
0.45 EE 0.99 0.10 E 0.78 0.16 FE 0.59 0.05 E 0.78 0.00 E 1.00 


(k) 土壤 湿度 指数 (2000 年 ) (1) 土壤 湿度 指数 (2005 年 ) (m) 土壤 湿度 指数 (2010 年 ) (n) 土壤 湿度 指数 (2015 年 ) (0) 土壤 湿度 指数 (2020 年 ) 


0 50 1 
土壤 湿度 指数 土壤 湿度 指数 土壤 湿度 指数 土壤 湿度 指数 土壤 湿度 指数 
0.00 0.93 0.09 f = 0.99 0.09 BM 0.95 0.00 F= 0.92 0.09 FE 0.99 


图 例 CO 黄河 源 园 区 边界 “三 3 水 域 及 其 边界 


图 6 土地 退化 影响 因子 分 布 
Fig. 6 Distributions of land degradation impact factors 


成 了 土壤 湿度 指数 与 沙化 特征 指数 对 第 下 类 的 空 
间 分 异 影 响 力 提 高 。 到 2020 年 ,降雨 量 在 偏 东 地 区 
有 所 增长 ,植被 覆盖 度 继 续 提 高 ,但 土壤 侵蚀 模 数 
整体 变化 不 大 。 随 着 气温 升 高 ,地 表 蒸 散 量 在 偏 北 
地 区 显著 提高 , 偏 北 地 区 的 土地 沙化 风险 增加 。 东 
北 侦 北 地 区 的 土壤 侵蚀 模 数 还 存在 小 幅 上 涨 ,造成 
了 土地 退化 程度 的 轻微 恶化 。 


4 讨论 


黄河 流域 生态 保护 与 高 质量 发 展 近 年 来 逐渐 
被 国家 重视 ,黄河 源头 更 应 得 到 更 多 的 关注 与 重 
视 。 研 究 结 果 指 出 ,降雨 侵蚀 力 因子 分 布 在 200~ 
300 MJ.mm.hm2.ha ,与 章 文 波 等 加 针对 全 国 所 
作 的 降雨 侵蚀 力 分 布 一 致 ;土壤 可 蚀 性 因子 分 布 在 
0.01~0.04 thm .hhn .MT mm 之 间 , 与 张 科 利 
等 ”的 研究 结果 一 致 ;其 余 因 子 均 来 自 官方 发 布 的 
DEM 数据 与 遥感 数据 直接 计算 ,并 未 进行 二 次 处 
理 。 同 时 , 胡 光 印 等 "对 黄河 流域 沙漠 化 问题 分 析 
发 现 ,黄河 源 园区 沙漠 化 问题 从 偏 东 南 地 区 向 偏 西 


北 地 区 逐渐 加 重 ,与 本 文 研究 结果 相符 。 

相 比 仅 依靠 土地 退化 指数 "进行 土地 退化 时 
空 演化 分 析 , 本 文 从 GeoAI 角 度 人 手 , 引 入 了 STICC 
聚 类 算法 "”。 聚 类 后 的 空间 分 异性 更 显著 ,与 许 明 
杰 等 外 研究 一 致 。 同 时 ,本 文 分 析 得 土地 退化 程度 
整体 呈现 波动 性 下 降 , 与 Ye 等 "针对 黄河 源 园区 的 
分 析 相 符 。 但 局 部 地 区 仍然 存在 问题 ,中 部 偏 南 地 
区 与 东北 偏 北 地 区 的 高 水 土 流失 问题 ,是 黄河 源 园 
区 土地 退化 空间 分 异 的 主要 原因 。 在 国家 政策 扶 
持 下 ,水土 流 失 问题 逐渐 缓解 ,土壤 侵蚀 模 数 对 空 
间 分 异 的 影响 力也 不 断 降 低 。 不 过 ,中 部 偏 北 地 区 
的 土地 沙化 风险 近 10 a 来 有 所 增加 ,主要 是 由 于 地 
表 莹 散 量 的 增加 ,这 与 气候 变 暖 息息相关 50。 特 别 
是 东北 地 区 ,水 土 流 失 风险 与 土地 沙化 风险 均 存 在 
小 幅 提 高 ,需要 引起 重视 。 


5 结论 


本 文 根 据 三 江 源 黄河 源 园区 2000—2020 年 存 
在 的 土地 退化 问题 ,综合 应 用 土壤 侵蚀 模 数 .沙化 
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特征 指数 与 土壤 湿度 指数 探究 了 该 区 域 土地 退化 
的 时 空 演化 过 程 ,主要 结论 如 下 : 

(1) 第 I 类 主要 分 布 在 中 部 偏 南 与 东北 地 区 ， 
属于 高 值 聚 类 区 ,分 布 面 积 较 少 ,以 水 土 流失 问题 
为 主 ;第 开 类 主要 分 布 在 偶 北 地 区 ,为 中 值 聚 类 区 ， 
多 年 分 布 面积 占 半数 左右 ,主要 表现 为 土地 沙化 问 
题 ;西南 与 东南 地 区 主要 分 布 为 土地 退化 程度 较 轻 
的 第 亚 类 

(2) 土地 退化 指数 整体 处 于 波动 性 下 降 , 近 20 a 
基本 呈现 轻 度 - 微 重 - 中 度 - 微 度 - 微 度 分 布 。 土 地 
退化 程度 在 黄河 源 园 区 整体 好 转 ,在 偏 南 地 区 的 好 
转 趋势 更 明显 。 

(3) 黄河 源 网 区 土地 退化 空间 分 异 现象 主要 受 
但 随 着 水 土 流失 问题 逐渐 被 治 

理 ,其 对 土地 退化 空间 分 异 的 影响 力 不 断 降低 。 同 
时 ,由 于 偏 北 地 区 近 些 年 来 渐渐 显示 出 土地 沙化 的 
趋势 ,该 区 域内 的 沙化 特征 指数 与 土壤 湿度 指数 对 
土地 退化 空间 分 异 的 影响 力 有 轻微 提高 。 
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Spatial-temporal evolution of land degradation in the Yellow River Source Park 
in recent 20 years 


QIN Tong', LI Gongquan', FAN Jiachen’ 
(1. School of Geosciences, Yangtze University, Wuhan 430000, Hubei, China; 2. School of Geography and Environmental Sciences, 
Northwest Normal University, Lanzhou 730070, Gansu, China) 


Abstract: Because Sanjiangyuan’ s Yellow River Source Park flows into the Yellow River, its ecosystem is 
crucial to the entire Yellow River Basin of China. Of the land degradation problems observed in the study area 
from 2000 to 2020, this study focuses on the soil erosion modulus, sanding characteristic index and soil moisture 
index. To achieve a more in-depth analysis of spatial and temporal evolution of land degradation in the area, machine 
learning algorithms were used. The spatial Toeplitz inverse covariance-based clustering method was used to carry 
out the cluster analysis of three parameters as part of the spatio-temporal evolution analysis of land degradation in 
the Yellow River Source Park. The extreme gradient boosting method was used to analyze the factors influencing 
the spatio-temporal evolution. The results of the study showed that: (1) The most significant land degradation was 
observed in the central-southern and northeastern regions, followed by the northern region. (2) In the northeast, 
there are two categories of problems: soil erosion and land desertification. In the central-southern region, severe 
soil erosion is the main problem. In the northern region, land desertification is the main problem. (3) In the past 
two decades, the land degradation index has shown a fluctuating decline that land degradation overall tends to 
improve, and that the effects are more obvious in the remote southern region. (4) The problem of soil erosion is 
improving overall, and there is a risk of land sanding in the remote northern areas. This study is a critical tool for 
understanding the spatial and temporal trends in the distribution and development of land degradation in the 
Yellow River Source Park and, by extension, for the promotion of increased ecological protection in Yellow River 
Source Park. 

Key words: land degradation; remote sensing; spatial Toeplitz inverse covariance- based clustering; extreme 


gradient boosting 


